NASA'S SATELLITE ORBIT ANOMALY PROBLEM 
CAN BE SOLVED PRECISELY 
IN THE FRAME OF EINSTEIN'S SPECIAL THEORY OF 

RELATIVITY. 

ANOMALY CONFIRMS THAT GRAVITY FIELDS PROPAGATE 
WITH VELOCITY OF LIGHT AS EINSTEIN PREDICTED 
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Abstract. NASA's Jet Propulsion Laboratory put on You Tube a problem 
that has been baffling the scientists for sometime. It involves an unexpected 
force acting on the space probes. 

The author proves that NASA'S satellite orbit anomaly problem can be 
solved in the frame of Einstein's Special Theory of Relativity. The anom- 
aly confirms that gravity fields propagate with velocity of light as Einstein 
predicted. 

The proof is based on the authors discovery of the relativistic version of 
Newton's gravity field. The author provides formulas for relativistic equation 
of motion for a spacecraft in the joint gravitational field of the Earth and the 
Sun in a Lorentzian frame attached to the Earth. The formulas are suitable 
for digital computers and can be easily implemented. He also shows how to 
find solutions of the relativistic equations of motion for the spacecraft. 



1. NASA IS BAFFLED BY AN UNEXPECTED FORCE 
ACTING ON SPACE PROBE 

The following is a transcript from You Tube video apparently made by Mr. 
Anderson from NASA's Jet Propulsion Laboratory pQ: 

Mysteriously five space craft that flew past the earth have each displayed unex- 
pected anomalies in their motions. Pioneer anomaly has hints that unexpected 
forces may appear to act on spacecraft. The anomalies were seen with identical Pi- 
oneer 10 and 11 spacecraft. Both seem to experience tiny but unexplained constant 
acceleration toward the sun. 

In five of the six flybys, the scientists have confirmed anomalies. "I'm feeling 
both humbled and perplexed" said Anderson. "There is something strange going 
on with spacecraft, with respect to Earth's equator. It suggests that the anomaly is 
related to Earth's rotation" 

"Another thing in common between the Pioneer and these flybys is what you 
would call an unbound orbit around a central body, " Anderson said. "For 
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instance, the Pioneers are flying out of the system. There 's definitely something 
going on. " 

"Could Earth's rotation be a part of the equation? Rotation forecasts are impor- 
tant to NASA 's Jet Propulsion Laboratory. " 

"The more forceful winds double the angular momentum of the atmosphere. 
Winds can effect the Earth 's rotation. " 

"Something's effecting the satellites orbit. Could it be Planet X? What do you 
think?" 

The purpose of this note is to answer the questions rased in this video and 
to provide the gravity field formulas in the frame of Einstein's special theory of 
relativity. These formulas can be used for navigation in gravitational fields in the 
presence of several celestial bodies in motion. 



2. The solution of the problem 

The answer comes from Einstein's special theory of relativity involving Maxwell's 
equations and formulas discovered by the author for relativistic gravity field that 
correct the Newton's formula. The anomaly indeed is due to Earth's rotation. 

With respect to the Lorentzian frame attached to Earth the Sun moves around 
the Earth almost periodically with period of 24 hours. The distance to the Sun is 
approximately 

r — 8 [mill] x c [km/ sec] 

so the circumference of Sun's orbit around the Earth in the Lorentzian frame is 
27rr cos(a), where a is the angle between the plane of Earth's equator and the 
plane of Sun's orbit, and the speed of the Sun in that frame is approximately 
4j cos(a) c. 
This fraction 

q = 30 C ° S ^ 

represents an important factor. It determines the rate of convergence in the com- 
putation of the field representing the retarded time. 

For the sake of simplicity select the unit of time so that the speed of light in 
the Lorentzian frame is c = 1. Select also the unit of mass so that gravitational 
constant G = 1. 

This means that if the unit of length is 1 meter, then the unit of time is approx- 
imately 3.3 nanoseconds and the unit of rest mass is about 3.8 metric tons. 

Assume that we have the formula t i— ► r(t) £ R 3 for the trajectory of the Sun 
around the Earth in our Lorentzian frame. It is almost periodic. So we may a priori 
estimate a bound on its velocity w(t) — r(t) in that frame. In our case it is q. 

Let the pair (y, t) £ R 3 x R denote any point in the Lorentzian frame. According 
to Einstein's theory of relativity [9] and his work jointly with Rosen 10J on gravity, 
gravity waves propagate through space with velocity of light. 

For a gravity wave from the Sun to arrive at a point y S R 3 at time t, the wave 
should be emitted from position r(r) on its trajectory at some earlier time t. This 
leads to the Lorentz relation 



t-r = \y- r(r)| 
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Since the Sun's trajectory is fixed the above equality implicitly defines the function 
(y, t) h-*- r. Notice that for given (y, t) the value r is a fixed point of the function 
f(s) =t — \y — r(s)\ considered on the interval s <t. The function / is a contraction 
and its Lipschitz constant is q. From Banach's fixed-point theorem, we get a fast 
converging algorithm for computing the function r(y,t). Namely use the recursive 
formula 

r o (y,t)=0 and r n (y, t) = /(r n _i(y, t) for all (y, t), n = 1, 2, . . . 
After n steps we get the following approximation 

n 

\T(y,t)-T n (y,t)\<f—\t-\y-r(0)\\ for all (y, t) e R\ n = 1, 2, . . . 

The function r is continuous as follows from results of Bogdan ,5], Theorem 3.1. 
Introduce the scalar field T called the time delay field by the formula 

T(y,t)=t-T(y,t) for all (y, t) G i? 4 

and vector field 

ri 2 (y,t) = y-r(T(y,t)) for all (y, t) e i? 4 

representing a vector starting on the trajectory of the Sun at retarded time and 
having its end at the point y £ R 3 at time t. It replaces Newton's radius vector. 
Notice the equality 

T(y,t) = \r 12 (y,t)\ 

for all (y,t) € i? 4 . 

Introduce also the velocity v and acceleration a vector fields by the formulas 

v(y,t) = w(r(y,t)) 
a(y,t) = w(T(y,t)) 

for all (y,t) € i? 4 . 
On the set 

G= {(y,t) £ R 4 : T(y,t)>0} 

introduce the scalar field u = 1/T, the unit vector field e = ur±2, representing 
direction of vector ri2, and the scalar field 

1 

z = . 

1 — e • v 

Notice that since the velocity field |v|<<7<c=lwe have for the dot product 

|e ■ v\ < q 

so the field z is well defined on the set G and we have the uniform estimate 

< z < — !— . 

" "1-9 

Definition 2.1 (Fundamental fields). The fields 

r, ri2, v, a, T, u, z, e 

will be called the fundamental fields associated with the trajectory of the moving 
point mass. 
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The fundamental fields arc continuous on their respective domains. This follows 
from the fact that composition of continuous functions yields a continuous function. 
Thus all of them, for sure, are continuous on the open set G of points that do not 
lie on the trajectory. 

We would like to stress here that the fundamental fields depend on the Lorentzian 
frame, in which we consider the trajectory. It is important to find expressions in- 
volving fundamental fields that yield fields invariant under Lorentzian transforma- 
tions. 

Lorentz and, independently, Einstein [5], Part II, section 6, established that 
fields satisfying Maxwell equations are invariant under Lorentzian transformations. 



Now consider the fields E, B defined by means of fundamental fields by the 
formulas 

E = u 2 e + u^D^e) + D 2 {e) 
B = e x E 

where D denotes the partial derivative with respect to time Jj. 

The pair (E, B) of fields associated with a trajectory of a moving point mass in a 
given Lorentzian frame satisfies homogenous Maxwell equations. Therefore accord- 
ing to a result of Lorentz and Einstein, the pair of fields (E, B) is invariant under 
Lorentzian transformations. The fields E and B, obtained purely by mathematical 
construction, form parts of an antisymmetric tensor of second rank. 

Indeed, the pair of fields 

E= (E 1 ,E 2 ,E 3 ) and B=(B 1 ,B 2 ,B 3 ) 

is a part of an antisymmetric tensor that in Lorentzian space-time i? 3 x R frame 
has the matrix that looks as follows 






+Ei 


+E 2 


+E 3 


Ex 





+B 3 


-B 2 


E 2 


-B 3 





+B X 


E 3 


+B 2 


-B x 






When the point mass is in its rest frame the field E differs from Newton's gravity 
field just by a constant of proportionality. The formula for the field E amends the 
formula derived heuristically by Feynman for a moving point charge. So it is proper 
to call the field E the Newton-Feynman field generated by the trajectory. 

The following theorem permits us to analyze such fields E and B by reducing 
the computations to simple algebraic expressions involving the fundamental fields. 

Introduce operators D = and Di = for i = 1, 2, 3 and V = (Dx, D 2 , D 3 ). 
Observe that Si in the following formulas denotes the i-th unit vector of the standard 
base in R 3 that is dx = (1,0,0), S 2 = (0,1,0), 5 3 = (0,0,1). 

Theorem 2.2 (Partial derivatives of fundamental fields). Assume that in some 
Lorentzian frame we are given an admissible trajectory t > r 2 (t). Define the time 
derivative w(t) = r 2 (t). For partial derivatives with respect to coordinates of the 
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vector n = y £ R 3 we have the following identities on the set G 



(2.1) 


D t T 




zei where v = ri o r, 




(2.2) 


DiU 




— zu 2 ei, 




(2.3) 


D iV 




—eiZa where a = w o r, 




(2.4) 


D,T 








(2.5) 


D ie 




—uzae + uSi + uzeiV where Si = 


(<5tf). 


(2.6) 


D lZ 




—z 3 ei(e, a) — uz 3 ei + uz 2 ei + uz 2 Vi 


f uz 3 


and for 


the partial derivative with respect to time we have 




(2.7) 




DT 


= 1-z, 




(2.8) 




Du 


2 2 

= zu — U , 




(2.9) 




Dt 


= z, 




(2.10) 




Dv 


= za where a = w o t, 




(2.11) 




De 


= —ue + uze — uzv, 




(2.12) 




Dz 


= uz — 2uz 2 + z 3 (e, a) + uz 3 — uz 


3 (v,v 



Since the expression on the right side of each formula represents a continuous func- 
tion, the fundamental fields are at least of class C 1 on the set G. Moreover if the 
trajectory is of class C°° then also the fundamental fields are of class C°° on G. 

3. Bogdan-Feynman Theorem for a moving point mass 



The following theorem is just a consequence that the wave emitted from a tra- 
jectory of a point mass propagates in a Lorentzian frame with the speed of light. 
Please notice that the physical nature of the fields is completely irrelevant. 

The important fact is that we are working in a Lorentzian frame, in which we 
are given a trajectory, on which the velocity and acceleration are bounded. Such 
trajectory alone generates in a unique way the system of the fundamental fields by 
means of which we are able to construct fields that are preserved by Lorentzian 
transformations. 

We remind the reader that we are working in a fixed Lorentzian frame. The 
position vector is denoted by r\ S R 3 and time by t € R. The position axes are 
oriented so as to form right screw orientation. The units are selected so that the 
speed of light c = 1. 

Partial derivatives with respect to coordinates of r\ are denoted by D\, D2, D3 
and with respect to time just by D. The gradient differential operator is denoted 
by V = (D 1 ,D 2 ,D 3 ) and the D'Alembertian operator by D 2 = V 2 - D 2 . The 
expression (e, v) denotes the dot product in R 3 of the vectors involved. 

Theorem 3.1 (Bogdan-Feynman Theorem). Assume that in a given Lorentzian 
frame the map t 1— * r 2 (t) from R to R 3 represents an admissible trajectory of class 
C 3 . Assume that G denotes the open set of points that do not lie on the trajectory. 
All the following field equations are satisfied on the entire set G. 

Consider the pair of fields E and B over the set G given by the formulas 

E = u 2 e + u~ 1 D(u 2 e) + D 2 e and B — e x E 

where u and e represent fundamental fields k2.1\) associated with the trajectory riii) . 
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Then this pair of fields will satisfy the following homogenous system of Maxwell 
equations 

VxE = -DB, V-E = 0, 

V x B = +DE, V • B = 0, 

and the homogenous wave equations 

U 2 E = 0, U 2 B = 0. 

Moreover Lienard- Wiechert potentials, expressed in terms of the fundamental 
fields as A = uzv and <f> — uz, satisfy the homogenous system of wave equations 
with Lorentz gauge formula 

n 2 A = o, n 2 <f> = o, v • a + D(p = o 

and generate the fields E and B by the formulas 

E = -V0 -DA and B = V x A. 

Finally we have the following explicit formula for the field E in terms of the 
fundamental fields 

E = —uz 2 a + uz 3 (e, a)e — uz 3 (e, a)v 

2 3 2 3/ \ 2S 23/ \ 

+ u z'' e — u z (v,v)e — u z v + u z (v,v)v. 

For a proof of the above theorem see Bogdan [5J, Theorem 11.1. As a consequence 
of the above theorem the components of the quantities E, B, A, and </>, since they 
satisfy the homogenous wave equation \3 2 y — 0, propagate in the Lorentzian frame 
with velocity of light c = 1. 

4. Relativistic equation of motion for a spacecraft 
IN the joint gravitational field 
of the Earth and the Sun 
in a Lorentzian frame attached to the Earth 

Now for a body of rest mass mo moving under the influence of a force F from 
Newton-Einstein formula, time derivative of the momentum is equal to force, we 
must have 

p = F 

where p = mo7 v is the relativistic momentum and 7(f) = (1— |u| 2 ) -1 ^ 2 . Computing 
the derivative of p with respect to time we get 

p = mo 77) + mo 7 v = mo 7 w + too 7 3 {v, v) v 

= m 7 (v + 7 2 (v, v) v) — too j(v) T(v)(v), 

where T(v) is the linear transformation of R 3 into R 3 given by the formula 

(4.1) I>)(/i) = h + j 2 (v,h)v for all h £ R 3 . 

Notice that for fixed velocity v the transformation T(v) represents a symmetric, 
positive definite transformation with eigenvalues equal respectively to 

(i + 7 2 M 2 ), 1, 1. 
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Thus the inverse transformation f (v) exists and is also symmetric. Since its eigen- 
values are 

(l+7 2 M 2 )-\ 1, 1 

and the norm of positive definite symmetric transformation is its maximal eigen- 
value, we must have 

|f»| = l for all veR 3 

and we have also 

t(v)T(v) = e for all v G R 3 

where e denotes here the identity transformation in the algebra Lin(R 3 , R 3 ) that 
is e(h) = h for all h G R 3 . 

We shall call the function v T(v) the reciprocal to the function T(v), since 

I» = [r(u)]- 1 for all v G R 3 ■ 

5. The mystery force identified 

Thus the relativistic equation of motion of the spacecraft can be written as 
(5.1) m j{v)T(v)(v) = -m e E e -m s E s 

where m e and m s denote the rest mass of the Earth and of the Sun, and where E e 
and E s represent the Newton-Fcynman fields generated by the Earth and the Sun 
respectively. The mystery force is F s = —m s E s and as Mr. Anderson suspected 
it is due to Earth's rotation, or equivalently due to the fact that in a Lorentzian 
frame attached to Earth the Sun moves around the Earth. 

6. Algorithm for computing the transformation Gamma hat 

The space Lin(U, U) of linear bounded operators from a Banach space U into 
itself beside being closed under addition and scalar multiplication is also closed 
under the composition of operators: P o Q G Lin(U, U) for all P, Q G Lin(U, U). 
This operation has the property that |PoQ|<|P||Q|. 

Definition 6.1 (Banach Algebras). A Banach space U is called a Banach alge- 
bra, if it is equipped with a bilinear operation (u, w) — > uw, from the product U xU 
into U, that is associative: (uw) z — u(w z), and such that 

\uw\ < \u\ \w\ for all u, w G U. 

If in addition there is an element e in U such that eu = ue = u for all u G U, 
then such an algebra is called a Banach algebra with unit. In such a case an 
element u is called invertible if for some w G U, called the inverse of u, we have 

uw = w u = e. 

The unit element and the inverse elements are unique. We denote the inverse of u 
by 
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The following proposition is very useful in numerical computations of the inverse 
transformations. It represents a simple application of the Banach fixed point theo- 
rem and it is essential in the development of the theory presented in this note. It 
represents Pro. 6.25 of Bogdan [4]. 

Proposition 6.2 (Inverse (e — t>) _1 exist if \v\ < 1). Let U be a Banach algebra 
with unit. Then for every element v G U such that \v\ < 1 the inverse w — (e — v)^ 1 
exists. It is a fixed point of the operator f defined by the following formula 

f{u) = e + vu for all iief/. 

Moreover the sequence w n — f(w n -±) where Wq = is of the form 

w n = e + v + v 2 H + v n — v k for all n > 

and we have the following estimate for the distance of the fixed point w and the 
approximation w n 

\v\ n 

\1V — W n \ < j— r for all n > 0, 

1 - M 

and we also have an explicit formula for the inverse element w as the sum of an 
absolutely convergent series 

w = e + v + v 2 ^ = v n . 

7. Returning to our problem at hand 

Notice first of all that the space Lin(R 3 ,R 3 ) of all linear transformations from 
R 3 into itself can be identified with the space of 3 x 3 matrices. If e denotes the 
unit matrix let H {v ) denote the transformation defined by the formula 

H(v)(h) = -7(V) 2 (v, h) v for all h G R 3 . 

The norm of this transformation is j~t^ts ■ So if the velocity \v\ stays below \/0.5 = 
0.707 . . . that is below the 70% of the speed of light the norm of the operator H 
will be less then 1. So we can use the above proposition to compute the inverse 
transformation 

f(v) = (e-H(v))- 1 

with any level of precision. 
8. Algorithm for computing relativistic trajectory of spacecraft 

So here we are. We can rewrite the equation of motion of the spacecraft (|5 . 1 [) in 
the equivalent form as 

V = {m -i(v)Y l T{v)(-m e E e -m s E s ) = f(y,t,v), 

V = v, 

and this system of differential equations we can solve, say, by means of Runge-Kutta 
methods with any level of precision. 

In closing I would like to thank Mr. Anderson from Jet Propulsion Laboratory 
for putting this problem onto You Tube. 
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Taking this opportunity I would like to thank also to my former colleagues, 
Donald Jezewski and Victor Bond, from Mission Planning and Analysis Division of 
NASA's Lyndon B. Johnson Space Center in Houston, Texas, for introducing me 
to problems of celestial mechanics. 
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